Anthropogenic land-use legacies underpin climate change-related risks to forest ecosystems.

Forest ecosystems with long-lasting human imprints can emerge worldwide as outcomes of land-use cessation. However, the interaction of these anthropogenic legacies with climate change impacts on forests is not well understood. Here, we set out how anthropogenic land-use legacies that persist in forest properties, following alterations in forest distribution, structure, and composition, can interact with climate change stressors. We propose a risk-based framework to identify anthropogenic legacies of land uses in forest ecosystems and quantify the impact of their interaction with climate-related stress on forest responses. Considering anthropogenic land-use legacies alongside environmental drivers of forest ecosystem dynamics will improve our predictive capacity of climate-related risks to forests and our ability to promote ecosystem resilience to climate change.

Functional traits and climate drive interspecific differences in disturbance-induced tree mortality.

With climate change, natural disturbances such as storm or fire are reshuffled, inducing pervasive shifts in forest dynamics. To predict how it will impact forest structure and composition, it is crucial to understand how tree species differ in their sensitivity to disturbances. In this study, we investigated how functional traits and species mean climate affect their sensitivity to disturbances while controlling for tree size and stand structure. With data on 130,594 trees located on 7617 plots that were disturbed by storm, fire, snow, biotic or other disturbances from the French, Spanish, and Finnish National Forest Inventory, we modeled annual mortality probability for 40 European tree species as a function of tree size, dominance status, disturbance type, and intensity. We tested the correlation of our estimated species probability of disturbance mortality with their traits and their mean climate niches. We found that different trait combinations controlled species sensitivity to disturbances. Storm-sensitive species had a high height-dbh ratio, low wood density and high maximum growth, while fire-sensitive species had low bark thickness and high P50. Species from warmer and drier climates, where fires are more frequent, were more resistant to fire. The ranking in disturbance sensitivity between species was overall consistent across disturbance types. Productive conifer species were the most disturbance sensitive, while Mediterranean oaks were the least disturbance sensitive. Our study identified key relations between species functional traits and disturbance sensitivity, that allows more reliable predictions of how changing climate and disturbance regimes will impact future forest structure and species composition at large spatial scales.

Global patterns of tree density are contingent upon local determinants in the world's natural forests.

Previous attempts to quantify tree abundance at global scale have largely neglected the role of local competition in modulating the influence of climate and soils on tree density. Here, we evaluated whether mean tree size in the world's natural forests alters the effect of global productivity on tree density. In doing so, we gathered a vast set of forest inventories including >3000 sampling plots from 23 well-conserved areas worldwide to encompass (as much as possible) the main forest biomes on Earth. We evidence that latitudinal productivity patterns of tree density become evident as large trees become dominant. Global estimates of tree abundance should, therefore, consider dependencies of latitudinal sources of variability on local biotic influences to avoid underestimating the number of trees on Earth and to properly evaluate the functional and social consequences.

Tree growth response to drought partially explains regional-scale growth and mortality patterns in Iberian forests.

Tree-ring data has been widely used to inform about tree growth responses to drought at the individual scale, but less is known about how tree growth sensitivity to drought scales up driving changes in forest dynamics. Here, we related tree-ring growth chronologies and stand-level forest changes in basal area from two independent data sets to test if tree-ring responses to drought match stand forest dynamics (stand basal area growth, ingrowth, and mortality). We assessed if tree growth and changes in forest basal area covary as a function of spatial scale and tree taxa (gymnosperm or angiosperm). To this end, we compared a tree-ring network with stand data from the Spanish National Forest Inventory. We focused on the cumulative impact of drought on tree growth and demography in the period 1981-2005. Drought years were identified by the Standardized Precipitation Evapotranspiration Index, and their impacts on tree growth by quantifying tree-ring width reductions. We hypothesized that forests with greater drought impacts on tree growth will also show reduced stand basal area growth and ingrowth and enhanced mortality. This is expected to occur in forests dominated by gymnosperms on drought-prone regions. Cumulative growth reductions during dry years were higher in forests dominated by gymnosperms and presented a greater magnitude and spatial autocorrelation than for angiosperms. Cumulative drought-induced tree growth reductions and changes in forest basal area were related, but initial stand density and basal area were the main factors driving changes in basal area. In drought-prone gymnosperm forests, we observed that sites with greater growth reductions had lower stand basal area growth and greater mortality. Consequently, stand basal area, forest growth, and ingrowth in regions with large drought impacts was significantly lower than in regions less impacted by drought. Tree growth sensitivity to drought can be used as a predictor of gymnosperm demographic rates in terms of stand basal area growth and ingrowth at regional scales, but further studies may try to disentangle how initial stand density modulates such relationships. Drought-induced growth reductions and their cumulative impacts have strong potential to be used as early-warning indicators of regional forest vulnerability.

Disentangling the Legacies of Climate and Management on Tree Growth.

Legacies of past climate conditions and historical management govern forest productivity and tree growth. Understanding how these processes interact and the timescales over which they influence tree growth is critical to assess forest vulnerability to climate change. Yet, few studies address this issue, likely because integrated long-term records of both growth and forest management are uncommon. We applied the stochastic antecedent modelling (SAM) framework to annual tree-ring widths from mixed forests to recover the ecological memory of tree growth. We quantified the effects of antecedent temperature and precipitation up to 4 years preceding the year of ring formation and integrated management effects with records of harvesting intensity from historical forest management archives. The SAM approach uncovered important time periods most influential to growth, typically the warmer and drier months or seasons, but variation among species and sites emerged. Silver fir responded primarily to past climate conditions (25-50 months prior to the year of ring formation), while European beech and Scots pine responded mostly to climate conditions during the year of ring formation and the previous year, although these responses varied among sites. Past management and climate interacted in such a way that harvesting promoted growth in young silver fir under wet and warm conditions and in old European beech under drier and cooler conditions. Our study shows that the ecological memory associated with climate legacies and historical forest management is species-specific and context-dependent, suggesting that both aspects are needed to properly evaluate forest functioning under climate change. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10021-021-00650-8.

Using spectral indices as early warning signals of forest dieback: The case of drought-prone Pinus pinaster forests.

Forest dieback processes linked to drought are expected to increase due to climate warming. Remotely sensed data offer several advantages over common field monitoring methods such as the ability to observe large areas on a systematic basis and monitoring their changes, making them increasingly used to assess changes in forest health. Here we aim to use a combined approximation of fieldwork and remote sensing to explore possible links between forest dieback and land surface phenological and trend variables derived from long Landsat time series. Forest dieback was evaluated in the field over 31 plots in a Mediterranean, xeric Pinus pinaster forest. Landsat 31-year time series of three greenness (EVI, NDVI, SAVI) and two wetness spectral indices (NMDI and TCW) were derived covering the period 1990-2020. Spectral indices from time series were decomposed into trend and seasonality using a Bayesian estimator while the relationships of the phenological and trend variables among levels of damage were assessed using linear and additive mixed models. We have not found any statistical pieces of evidence of extension or shortening patterns for the length of the phenological season over the examined 31-year period. Our results indicate that the dieback process was mainly related to the trend component of the spectral indices series whereas the phenological metrics were not related to forest dieback. We also found that plots with more dying or damaged trees displayed lower spectral indices trends after a severe drought event in the middle of the 1990s, which confirms the Landsat-derived spectral indices as indicators of early-warning signals. Drops in trends occurred earlier for wetness indices rather than for greenness indices which suggests that the former could be more appropriate for dieback detection, i.e. they could be used as early warning signals of impending loss of tree vigor.

Climate reverses directionality in the richness-abundance relationship across the World's main forest biomes.

More tree species can increase the carbon storage capacity of forests (here referred to as the more species hypothesis) through increased tree productivity and tree abundance resulting from complementarity, but they can also be the consequence of increased tree abundance through increased available energy (more individuals hypothesis). To test these two contrasting hypotheses, we analyse the most plausible pathways in the richness-abundance relationship and its stability along global climatic gradients. We show that positive effect of species richness on tree abundance only prevails in eight of the twenty-three forest regions considered in this study. In the other forest regions, any benefit from having more species is just as likely (9 regions) or even less likely (6 regions) than the effects of having more individuals. We demonstrate that diversity effects prevail in the most productive environments, and abundance effects become dominant towards the most limiting conditions. These findings can contribute to refining cost-effective mitigation strategies based on fostering carbon storage through increased tree diversity. Specifically, in less productive environments, mitigation measures should promote abundance of locally adapted and stress tolerant tree species instead of increasing species richness.

Evidence of non-stationary relationships between climate and forest responses: Increased sensitivity to climate change in Iberian forests.

Climate and forest structure are considered major drivers of forest demography and productivity. However, recent evidence suggests that the relationships between climate and tree growth are generally non-stationary (i.e. non-time stable), and it remains uncertain whether the relationships between climate, forest structure, demography and productivity are stationary or are being altered by recent climatic and structural changes. Here we analysed three surveys from the Spanish Forest Inventory covering c. 30 years of information and we applied mixed and structural equation models to assess temporal trends in forest structure (stand density, basal area, tree size and tree size inequality), forest demography (ingrowth, growth and mortality) and above-ground forest productivity. We also quantified whether the interactive effects of climate and forest structure on forest demography and above-ground forest productivity were stationary over two consecutive time periods. Since the 1980s, density, basal area and tree size increased in Iberian forests, and tree size inequality decreased. In addition, we observed reductions in ingrowth and growth, and increases in mortality. Initial forest structure and water availability mainly modulated the temporal trends in forest structure and demography. The magnitude and direction of the interactive effects of climate and forest structure on forest demography changed over the two time periods analysed indicating non-stationary relationships between climate, forest structure and demography. Above-ground forest productivity increased due to a positive balance between ingrowth, growth and mortality. Despite increasing productivity over time, we observed an aggravation of the negative effects of climate change and increased competition on forest demography, reducing ingrowth and growth, and increasing mortality. Interestingly, our results suggest that the negative effects of climate change on forest demography could be ameliorated through forest management, which has profound implications for forest adaptation to climate change.

How do trees respond to species mixing in experimental compared to observational studies?

For decades, ecologists have investigated the effects of tree species diversity on tree productivity at different scales and with different approaches ranging from observational to experimental study designs. Using data from five European national forest inventories (16,773 plots), six tree species diversity experiments (584 plots), and six networks of comparative plots (169 plots), we tested whether tree species growth responses to species mixing are consistent and therefore transferrable between those different research approaches. Our results confirm the general positive effect of tree species mixing on species growth (16% on average) but we found no consistency in species-specific responses to mixing between any of the three approaches, even after restricting comparisons to only those plots that shared similar mixtures compositions and forest types. These findings highlight the necessity to consider results from different research approaches when selecting species mixtures that should maximize positive forest biodiversity and functioning relationships.

Disentangling the relative role of climate change on tree growth in an extreme Mediterranean environment.

Climate change can impair ecosystem functions and services in extensive dry forests worldwide. However, attribution of climate change impacts on tree growth and forest productivity is challenging due to multiple inter-annual patterns of climatic variability associated with atmospheric and oceanic circulations. Moreover, growth responses to rising atmospheric CO2, namely carbon fertilization, as well as size ontogenetic changes can obscure the climate change signature as well. Here we apply Structural Equation Models (SEM) to investigate the relative role of climate change on tree growth in an extreme Mediterranean environment (i.e., extreme in terms of the combination of sandy-unconsolidated soils and climatic aridity). Specifically, we analyzed potential direct and indirect pathways by which different sources of climatic variability (i.e. warming and precipitation trends, the North Atlantic Oscillation, [NAO]; the Mediterranean Oscillation, [MOI]; the Atlantic Mediterranean Oscillation, [AMO]) affect aridity through their control on local climate (in terms of mean annual temperature and total annual precipitation), and subsequently tree productivity, in terms of basal area increments (BAI). Our results support the predominant role of Diameter at Breast Height (DHB) as the main growth driver. In terms of climate, NAO and AMO are the most important drivers of tree growth through their control of aridity (via effects of precipitation and temperature, respectively). Furthermore and contrary to current expectations, our findings also support a net positive role of climate warming on growth over the last 50 years and suggest that impacts of climate warming should be evaluated considering multi-annual and multi-decadal periods of local climate defined by atmospheric and oceanic circulation in the North Atlantic.

Erratum: Complementarity effects on tree growth are contingent on tree size and climatic conditions across Europe.

This corrects the article DOI: 10.1038/srep32233.

Distribution of pines in the Iberian Peninsula agrees with species differences in foliage frost tolerance, not with vulnerability to freezing-induced xylem embolism.

Drought and frosts are major determinants of plant functioning and distribution. Both stresses can cause xylem embolism and foliage damage. The objective of this study was to analyse if the distribution of six common pine species along latitudinal and altitudinal gradients in Europe is related to their interspecific differences in frost tolerance and to the physiological mechanisms underlying species-specific frost tolerance. We also evaluate if frost tolerance depends on plant water status. We studied survival to a range of freezing temperatures in 2-year-old plants and assessed the percentage loss of hydraulic conductivity (PLC) due xylem embolism formation and foliage damage determined by needle electrolyte leakage (EL) after a single frost cycle to -15 °C and over a range of predawn water potential (ψpd) values. Species experiencing cold winters in their range (Pinus nigra J.F. Arnold, Pinus sylvestris L. and Pinus uncinata Raymond ex A. DC.) had the highest frost survival rates and lowest needle EL and soluble sugar (SS) concentration. In contrast, the pines inhabiting mild or cool winter locations (especially Pinus halepensis Mill. and Pinus pinea L. and, to a lesser extent, Pinus pinaster Ait.) had the lowest frost survival and highest needle EL and SS values. Freezing-induced PLC was very low and differences among species were not related to frost damage. Reduction in ψpd decreased leaf frost damage in P. pinea and P. sylvestris, increased it in P. uncinata and had a neutral effect on the rest of the species. This study demonstrates that freezing temperatures are a major environmental driver for pine distribution and suggests that interspecific differences in leaf frost sensitivity rather than vulnerability to freezing-induced embolism or SS explain pine juvenile frost survival.

Last-century forest productivity in a managed dry-edge Scots pine population: the two sides of climate warming.

Climate change in the Mediterranean, associated with warmer temperatures and more frequent droughts, is expected to impact forest productivity and the functioning of forests ecosystems as carbon reservoirs in the region. Climate warming can positively affect forest growth by extending the growing season, whereas increasing summer drought generally reduces forest productivity and may cause growth decline, trigger dieback, hamper regeneration, and increase mortality. Forest management could potentially counteract such negative effects by reducing stand density and thereby competition for water. The effectiveness of such interventions, however, has so far mostly been evaluated for short time periods at the tree and stand levels, which limits our confidence regarding the efficacy of thinning interventions over longer time scales under the complex interplay between climate, stand structure, and forest management. In this study, we use a century-long historical data set to assess the effects of climate and management on forest productivity. We consider rear-edge Scots pine (Pinus sylvestris) populations covering continental and Mediterranean conditions along an altitudinal gradient in Central Spain. We use linear mixed-effects models to disentangle the effects of altitude, climate, and stand volume on forest growth and ingrowth (recruitment and young trees' growth). We find that warming tends to benefit these tree populations, warmer winter temperature has a significant positive effect on both forest growth and ingrowth, and the effect is more pronounced at low elevations. However, drought conditions severely reduce growth and ingrowth, in particular when competition (stand volume) is high. We conclude that summer droughts are the main threat to Scots pine populations in the region, and that a reduction of stand volume can partially mitigate the negative impacts of more arid conditions. Mitigation and adaptation measures could therefore manage stand structure to adopt for the anticipated impacts of climate change in Mediterranean forest ecosystems.

Continental mapping of forest ecosystem functions reveals a high but unrealised potential for forest multifunctionality.

Humans require multiple services from ecosystems, but it is largely unknown whether trade-offs between ecosystem functions prevent the realisation of high ecosystem multifunctionality across spatial scales. Here, we combined a comprehensive dataset (28 ecosystem functions measured on 209 forest plots) with a forest inventory dataset (105,316 plots) to extrapolate and map relationships between various ecosystem multifunctionality measures across Europe. These multifunctionality measures reflected different management objectives, related to timber production, climate regulation and biodiversity conservation/recreation. We found that trade-offs among them were rare across Europe, at both local and continental scales. This suggests a high potential for 'win-win' forest management strategies, where overall multifunctionality is maximised. However, across sites, multifunctionality was on average 45.8-49.8% below maximum levels and not necessarily highest in protected areas. Therefore, using one of the most comprehensive assessments so far, our study suggests a high but largely unrealised potential for management to promote multifunctional forests.

Forest productivity in southwestern Europe is controlled by coupled North Atlantic and Atlantic Multidecadal Oscillations.

The North Atlantic Oscillation (NAO) depicts annual and decadal oscillatory modes of variability responsible for dry spells over the European continent. The NAO therefore holds a great potential to evaluate the role, as carbon sinks, of water-limited forests under climate change. However, uncertainties related to inconsistent responses of long-term forest productivity to NAO have so far hampered firm conclusions on its impacts. We hypothesize that, in part, such inconsistencies might have their origin in periodical sea surface temperature anomalies in the Atlantic Ocean (i.e., Atlantic Multidecadal Oscillation, AMO). Here we show strong empirical evidence in support of this hypothesis using 120 years of periodical inventory data from Iberian pine forests. Our results point to AMO+ NAO+ and AMO-NAO- phases as being critical for forest productivity, likely due to decreased winter water balance and abnormally low winter temperatures, respectively. Our findings could be essential for the evaluation of ecosystem functioning vulnerabilities associated with increased climatic anomalies under unprecedented warming conditions in the Mediterranean.

Climate- and successional-related changes in functional composition of European forests are strongly driven by tree mortality.

Intense droughts combined with increased temperatures are one of the major threats to forest persistence in the 21st century. Despite the direct impact of climate change on forest growth and shifts in species abundance, the effect of altered demography on changes in the composition of functional traits is not well known. We sought to (1) quantify the recent changes in functional composition of European forests; (2) identify the relative importance of climate change, mean climate and forest development for changes in functional composition; and (3) analyse the roles of tree mortality and growth underlying any functional changes in different forest types. We quantified changes in functional composition from the 1980s to the 2000s across Europe by two dimensions of functional trait variation: the first dimension was mainly related to changes in leaf mass per area and wood density (partially related to the trait differences between angiosperms and gymnosperms), and the second dimension was related to changes in maximum tree height. Our results indicate that climate change and mean climatic effects strongly interacted with forest development and it was not possible to completely disentangle their effects. Where recent climate change was not too extreme, the patterns of functional change generally followed the expected patterns under secondary succession (e.g. towards late-successional short-statured hardwoods in Mediterranean forests and taller gymnosperms in boreal forests) and latitudinal gradients (e.g. larger proportion of gymnosperm-like strategies at low water availability in forests formerly dominated by broad-leaved deciduous species). Recent climate change generally favoured the dominance of angiosperm-like related traits under increased temperature and intense droughts. Our results show functional composition changes over relatively short time scales in European forests. These changes are largely determined by tree mortality, which should be further investigated and modelled to adequately predict the impacts of climate change on forest function.

Complementarity effects on tree growth are contingent on tree size and climatic conditions across Europe.

Neglecting tree size and stand structure dynamics might bias the interpretation of the diversity-productivity relationship in forests. Here we show evidence that complementarity is contingent on tree size across large-scale climatic gradients in Europe. We compiled growth data of the 14 most dominant tree species in 32,628 permanent plots covering boreal, temperate and Mediterranean forest biomes. Niche complementarity is expected to result in significant growth increments of trees surrounded by a larger proportion of functionally dissimilar neighbours. Functional dissimilarity at the tree level was assessed using four functional types: i.e. broad-leaved deciduous, broad-leaved evergreen, needle-leaved deciduous and needle-leaved evergreen. Using Linear Mixed Models we show that, complementarity effects depend on tree size along an energy availability gradient across Europe. Specifically: (i) complementarity effects at low and intermediate positions of the gradient (coldest-temperate areas) were stronger for small than for large trees; (ii) in contrast, at the upper end of the gradient (warmer regions), complementarity is more widespread in larger than smaller trees, which in turn showed negative growth responses to increased functional dissimilarity. Our findings suggest that the outcome of species mixing on stand productivity might critically depend on individual size distribution structure along gradients of environmental variation.

Jack-of-all-trades effects drive biodiversity-ecosystem multifunctionality relationships in European forests.

There is considerable evidence that biodiversity promotes multiple ecosystem functions (multifunctionality), thus ensuring the delivery of ecosystem services important for human well-being. However, the mechanisms underlying this relationship are poorly understood, especially in natural ecosystems. We develop a novel approach to partition biodiversity effects on multifunctionality into three mechanisms and apply this to European forest data. We show that throughout Europe, tree diversity is positively related with multifunctionality when moderate levels of functioning are required, but negatively when very high function levels are desired. For two well-known mechanisms, 'complementarity' and 'selection', we detect only minor effects on multifunctionality. Instead a third, so far overlooked mechanism, the 'jack-of-all-trades' effect, caused by the averaging of individual species effects on function, drives observed patterns. Simulations demonstrate that jack-of-all-trades effects occur whenever species effects on different functions are not perfectly correlated, meaning they may contribute to diversity-multifunctionality relationships in many of the world's ecosystems.

Biotic homogenization can decrease landscape-scale forest multifunctionality.

Many experiments have shown that local biodiversity loss impairs the ability of ecosystems to maintain multiple ecosystem functions at high levels (multifunctionality). In contrast, the role of biodiversity in driving ecosystem multifunctionality at landscape scales remains unresolved. We used a comprehensive pan-European dataset, including 16 ecosystem functions measured in 209 forest plots across six European countries, and performed simulations to investigate how local plot-scale richness of tree species (α-diversity) and their turnover between plots (ß-diversity) are related to landscape-scale multifunctionality. After accounting for variation in environmental conditions, we found that relationships between α-diversity and landscape-scale multifunctionality varied from positive to negative depending on the multifunctionality metric used. In contrast, when significant, relationships between ß-diversity and landscape-scale multifunctionality were always positive, because a high spatial turnover in species composition was closely related to a high spatial turnover in functions that were supported at high levels. Our findings have major implications for forest management and indicate that biotic homogenization can have previously unrecognized and negative consequences for large-scale ecosystem multifunctionality.